The Trojan Horse Method: an Indirect Technique in Nuclear Astrophysics. 
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The Trojan Horse (TH) method is a powerful indirect technique that provides information to 
determine astrophysical factors for rearrangement processes at astrophysically relevant energies. A 
short coming for understanding the reliability of the technique has been determining the importance 
of nuclear and Coulomb effects on the energy dependence of the yield. Using a simple model, we 
demonstrate that off-energy-shell and Coulomb effects in the entry channel and the final state nuclear 
interactions do not change the energy dependence of the astrophysical factor extracted from the TH 
reaction. Some examples are presented. 



PACS numbers: 26.20.+I, 24.50.+g, 25.70.Ef, 25.70.Hi 

The presence of the Coulomb barrier for colliding 
charged nuclei makes nuclear reaction cross sections at 
astrophysical energies so small that their direct measure- 
ment in the laboratory is very difficult, or even impos- 
sible. Consequently indirect techniques often are used 
to determine these cross sections. The Trojan Horse 
(TH) method is a powerful indirect technique which al- 
lows one to determine the astrophysical factor for rear- 
rangement reactions. The TH method, first suggested by 
Baur [|, involves obtaining the cross section of the bi- 
nary x + A — > b + B process at astrophysical energies by 
measuring the two-body to three-body (2 — > 3) process, 
a + A — > y + b + B, in the quasifree (QF) kinematics 
regime, where the "Trojan Horse" particle, a = (xy), is 
accelerated at energies above the Coulomb barrier. Af- 
ter penetrating through the Coulomb barrier, nucleus a 
undergoes breakup leaving particle x to interact with tar- 
get A while projectile y flies away. From the measured 
a + A^y + b + B cross section, the energy dependence 
of the binary subprocess, x + A — > b + B, is determined. 

The main advantage of the TH method is that the ex- 
tracted cross section of the binary subprocess does not 
contain the Coulomb barrier factor. Consequently the 
TH cross section can be used to determine the energy de- 
pendence of the astrophysical factor, S(E), of the binary 
process, x + A — > b + B, down to zero relative kinetic 
energy of the particles x and A without distortion due 
to electron screening 2, 3]. The absolute value of S{E) 
must be found by normalization to direct measurements 
at higher energies. At low energies where electron screen- 
ing becomes important, comparison of the astrophysical 
factor determined from the TH method to the direct re- 
sult provides a determination of the screening potential. 

Even though the TH method has been applied success- 
fully to many direct and resonant processes (see Q and 
references therein), there are still reservations about the 
reliability of the method due to two potential modifica- 
tions of the yield from off-shell effects and initial and final 
state interactions in the TH 2^3 reaction. In the TH 
reaction, shown schematically in Fig. ^ particle x in the 
binary subprocess x + A — > b + B is virtual (off-energy- 
shell). In the standard analysis, the virtual nature of x is 




FIG. 1: The pole diagram describing the TH reaction a+A — > 
y + b + B. 

neglected and the plane wave approximation is used 
Here we address, for the first time, the reliability of this 
assumption. We also consider scattering between parti- 
cles a and A in the initial channel of the TH reaction and 
the dominance of the QF mechanism. Note that the TH 
reaction is a many-body process (at least four-body) and 
its strict analysis requires many-body techniques. How- 
ever some important features of the TH method can be 
addressed in a simple model. 

First the notation for the subsequent discussion is 
introduced. Eij{= kfj/2/j,ij) and are the relative 
energy and momentum of the real particles i and j, 
Hij = mi m,j/{mi + m,j), where rrii is the mass of particle 
i, and p x j is the relative momentum of the virtual par- 
ticle x and particle j. The internal particle x of the dia- 
gram shown in Fig. ^is virtual, i.e. E x ^ P x /2 m x . From 
energy-momentum conservation laws in the three-ray ver- 
tex a — ► x + y and the four-ray vertex x + A — > b + B, 
we get a x = p 2 xy /{2^ xy ) + e a xy = p 2 xA /(2fi xA ) - E xA , 
where e xy = m x + m y — m a is the binding energy for 
the virtual decay a — > x + y. Thus in a TH reaction 
p 2 xA /(2p xA ) > E xA = k 2 xA /{2p xA ) holds. The reaction 
amplitude corresponding to the diagram is given by (all 
particles are assumed, for simplicity, to be spinless) M p = 

M 2 ( 5? F) VxyiPxy)- Here M^° F) {a x , E xA , z) is the half- 
off-energy-shell (HOF) reaction amplitude for the direct 
binary subprocess x+A — ► b+B which depends on the ad- 
ditional variable a x due to the virtual nature of particle x. 
Also z = Pz^-kbs, with k = k/fc. The Fourier transform 
of the s-wave bound-state wave function for a — (xy) can 
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be written as tp xy (p xy ) = —W xy {p X y)/(^ x ) 1 nxy , where 
rj xy is the Coulomb parameter of the bound state (xy) 
and W xy is the amplitude for the virtual decay a — > x + y, 
which is regular at the singular point a x — 0. This sin- 
gularity is a branch point for the charged particles x and 
y. For a = d = (pn), (p X y(p xy ) has a pole at a x = 0, 
corresponding to the real (on-energy-shell (ON)) particle 
x. Hence the name pole diagram for Fig. ^ 6]. The 
modulus of the amplitude \M p \ of the pole diagram has 
a maximum at p xy — called the QF peak and the con- 
dition corresponding to p xy — is called QF kinematics. 
The QF peak is a trace of the pole located at p xy < 0. 
The QF kinematics are ideal for the TH method since 
small p xy corresponds to large separation distance be- 
tween particles x and y, thereby allowing particle y to be 
treated as a spectator. 

By dropping all the terms containing particle y and its 
interaction with other nuclei, the half-off-shcll post-form 
amplitude of the direct binary reaction, A(x, b)B, can be 
extracted from the exact amplitude for the TH reaction, 
A(a, yb)B. The result is given by 
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(HOF), 



(a x ,E xA ,z) ~< Pb<PB\d\<p x (pAP x A > 



(1) 

Here, ip^g is the distorted wave in the final state of the 
binary process, <pi is the bound state wave function of 
nucleus i, O — AVbB [1 + G^2 AV x a] is the transition 
operator, AVij — Vij — Uij, Vij {Uij) is the interaction 

(optical) potential between nuclei i and j and G^2 IS 
the Green's function of the system x + A (or b + B). The 
half-off-shell amplitude contains the off-shell plane wave 
\p x A > which describes the relative motion of the vir- 
tual particle x and A in the initial channel of the binary 
reaction rather than the distorted wave describing the 
initial state in the on-shell amplitude. Hence the half- 
off-shell amplitude does not contain a Coulomb barrier 
factor. First we estimate the contribution from the tran- 
sition operator [A VbB G^2 A V x a] to M^^ F ' . Consider 
two reactions previously analyzed using the TH method, 
6 Li(d,a) 4 He (Q 2 ^2 = 21.64 MeV) and 7 Li(p,a) 4 He 
(Q2^2 = 16.56 MeV) Q- The first reaction was treated 
as a deuteron transfer and the second as a triton transfer. 
Nuclei A and b are bound states A = (Bt) and b = (xt) 
with constituent particles x, t, B, where t is the trans- 
ferred particle. In this simple model the half-off-shcll 
amplitude takes the form (all particles are assumed to be 
spinless) 



M Z^2 (<T X ,E x a;z) ~< tpl B ' (p x t\0\(pBtP X A > ■ (2) 

The bound state wave functions have been approximated 
by their tails, AVbB ~ V B t and the zero-range approxi- 
mation (ZRA) has been used for Vst^Bt- The ZRA is 
good enough to determine the energy dependence of the 
cross section. At E x a < 1 MeV, only the Coulomb in- 
teraction needs to be included in the transition operator 
AV x a and in the Green's function. Also at low E xA the 
dominant contribution comes only from the s-wave in the 
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FIG. 2: (Color online). Energy dependence (E = E x a) of 
the half-off-shell (red dashed line) and on-shell (black solid 
line) astrophysical factors for (a) the 7 Li(p, a) 4 He reaction 
(a = d, x — p, y = n); (b) the 6 Li(d, a) 4 He reaction (a = 
6 Li, x — d, y — a). 



channel x + A due to the large Q^^i value for both re- 
actions. Consequently the z dependence of M^^^ can 
be neglected. The contribution from the transition op- 
erator [AVbBG x + 2 AV x a] is negligible for both reactions 
as is the case for the on-shell processes. It constitutes a 
background for the DWBA with the transition operator 
given by AV b B- 

The on-shell amplitude to be compared with the half- 
off-shell one is 



MiZ N 2 \E xA ) ~< ^B^AVbB^Bt^ 



> 



(3) 



Here, V^aj ^ ne distorted wave describing the scatter- 
ing of real particles x and A in the initial channel, 
can be approximated by the pure Coulomb scattering 
wave function at low energies. Consequently, ip x + 2 — 
N x a exp(i \i x A • ? x a) iFi, where N x a is the Gamow nor- 
malization factor and \F\ is the hypergeometric function, 
which has a very weak energy dependence at small k X A- 
The only difference in the energy dependence between 
the half- off- shell and on-shell astrophysical factors comes 
from the the use of the half-off-shcll plane wave in the 
initial state in Eq. J2J and the on-shell distorted wave in 
Eq. From Fig. [21 it is clear that the results, which 
were calculated for QF kinematics, clearly justify the TH 
method. The energy dependence of the half-off-shell and 
on-shell astrophysical factors at low energies are practi- 
cally identical. Since only the energy dependence is of 
interest, the half-off-shell result in Fig. [21 has been nor- 
malized to the on-shell one at an energy E x a = 1 keV 
for ease of comparison. 

Another way to determine the effect of the virtual na- 
ture of the transferred particle is to vary the parameter 
a x . In Fig. [21 (a) the energy dependence of the half-off- 
shell astrophysical factor for the 6 Li(d, a) 4 He reaction 
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FIG. 3: (Color online), (a) Energy dependence (E = E d a Li ) 
of the half-off-shell astrophysical factor for the 6 Li(d, a) 4 He 
reaction as a function of as (?d — 0- black solid line (on- 
shell kinematics); ad = £da~ re d dashed line (QF kinematics); 
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1-5 £ da ~ Dme dotted line; ad = 2.5 e da - green dashed- 



dotted line, (b) Comparison of the energy dependence of the 
half-off-shell astrophysical factor for the 6 Li(d, a) 4 He reaction 
determined from the 6 Li( 6 Li, a a) 4 He TH reaction with (black 
solid line) and without (red dashed line) inclusion of the initial 
6 Li — 6 Li Coulomb interaction, (c) The same as in (a) but 
with the final-state interaction described by the Ali-Bodmer 
potential: ad = 0- black solid line (on-shell kinematics) ; ad = 
£ <ia' re d dashed line (QF kinematics). In both cases (b) and 
(c) the red dashed line is normalized to the black solid line at 
E = 1 keV. 

(x = d; A = 6 Li; b = B = a) is plotted for different 
values of ad- Note that — corresponds to the on- 
shell deuteron. The deuteron becomes farther from the 
energy-shell as grows. It is clear from Fig. [3] that 
the size of the S factor changes but the energy depen- 
dence does not change. This justifies the procedure of 
disregarding the virtual nature of the entryparticle in 
the analysis of the binary process A(x,b)B If the 

amplitude of the TH process A(a, y b)B is approximated 
by the pole amplitude M p , as is done in the TH method, 
then the amplitude of the binary subprocess can be triv- 
ially singled out. This approximation ignores the inter- 
action between particles a and A in the entry channel. 
Inclusion of the Coulomb interaction in the a — A chan- 
nel is a way to gauge the validity of this approximation. 
Note that including the Coulomb-nuclear interaction in 
the exit channel does not change the final result. The 
TH reaction amplitude with Coulomb scattering in the 
entrance channel is given by 



TH) 



(4) 



where ip y F 1S ^ ne distorted wave describing the scattering 
of y in the center-of-mass of the system F = b + B. As an 
example, consider the TH reaction 7 Li(d, n cv) 4 He, where 
a = d = (pn), A = 7 Li. The bound state wave function 



of d is approximated by its tail. The Coulomb interaction 
between d— 7 Li modifies the behavior of the amplitude of 
the TH reaction. For simplicity, the final state scattering 
wave functions are replaced by plane waves. To single 
out the amplitude of the binary process 7 Li(p, cv) 4 He in 
QF kinematics, the M^ TH ^ is divided by the factor 



R =< IpxyQtyF ~ (m y / m a ) p aA ) | Ip^J} (Ptul) >, 

performed over p a A 



(5) 



where integration is performed over p a A and 
(Po,a) is the Fourier component of the a — A 
Coulomb scattering wave function. This factor de- 
termines the behavior of M^ TH ^ near the singularity 
a x = 0. Note that for the results presented in Fig. [2] the 
half-off-shell astrophysical factor has been calculated at 
a fixed p xy — corresponding to QF kinematics of the 
TH 2 — > 3 reaction. E x a becomes a unique function 
of E a A- In practice TH experiments are carried out 
at a fixed energy, E a A, while p xy is allowed to vary 
around p xy = 0. For example, to cover the energy region 
< E p 7 U < 800 keV for the 7 Li(d, na) 4 He reaction at 
energy E d r u = 4.20 MeV, it is sufficient to vary the 
factor <7p by 7% around the QF value a p = 2.224 MeV. 

To test the importance of the Coulomb interac- 
tion in the initial state, the energy dependence of 
the half-off-shell astrophysical factor for the binary 
6 Li(d, a) 4 He reaction, determined from the TH reaction 
6 Li( 6 Li,aa) 4 He at & Li 6 Li = 3.14 MeV, the TH reac- 
tion amplitude (Eq. was calculated with and without 
the initial Coulomb interaction term. The results, as 
seen in Fig. [3] (b), show that the energy dependence of 
both astrophysical factors is nearly identical. Thus in- 
cluding the Coulomb interaction in the initial state does 
not affect the energy dependence of the astrophysical fac- 
tor determined from the TH reaction. Similar calcula- 
tions done for S l (£ , p 7 L i) for the binary 7 Li(p, a) 4 He reac- 
tion, determined from the TH reaction Li(d, n a) He at 
Ed 7 u — 4.03 MeV, show the same result. As also seen 
from Fig. 01(c), the final state interaction (here described 
by the Ali-Bodmer a — a potential) does not affect the 
energy behavior of the on-shell and half-off-shell astro- 
physical factors for the binary process 6 Li(d, a) 4 He. 

Understanding TH reactions proceeding through a res- 
onant binary subprocess is simpler than for the direct 
case. The TH pole mechanism for the direct process 
near the QF peak is dominated by the singularity (pole 
for a = d or branch point for charged particles x and y) 
in the bound state wave function <p xy {p xy )- In a reso- 
nant reaction, the mechanism depends on two singulari- 
ties, one that occurs in the bound state wave function in 
the p xy plane and the other that is due to the resonant 
pole of the amplitude of the binary subprocess A(x, b)B 
in the E x a plane. The TH 2^3 reaction proceeding 
through the resonance in the subsystem x + A is very 
similar to stripping (d,p), (d,n) or a-transfer ( 6 Li, d) re- 
actions populating resonant states The difference is 
that in the TH process we are interested in the decay of 
the intermediate resonances into the channel b + B which 
differs from the entry channel x + A. The half-off-shell 
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FIG. 4: S(£) for the 7 Li(p, a) 4 He reaction. The solid curve 
is the polynomial fit to the data (full circles) obtained from 
the TH 7 Li(d, an) 4 He reaction. Open cirles are the data 
from the direct measurement which include electron screening 
effects H. 



s-wave resonant amplitude is given by 



M (»OF )W{kb ^ pxA . E) = _ 27T 
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J xA 



(6) 



Here, TbB{EbB) is the partial width of the resonance 
F* in the channel b + B, k 



xA(R) 



^xA IS ^ ne resonance energy and Sfo is the non- 
resonant (potential) scattering phase shift of particles 
b and B in the final state. The off-shell s-wave 
form factor for the vertex x + A — > F* is given 

hy W (p xA ,k x A(B)) = [E { X A ~ P 2 xAl{ 2 ^A)]i>'^ ) (pxA), 

where V„o (Pxa) =< jo(PxAr)\^' (r) > is the Fourier 
transform of the Gamow radial wave function Vw) ( r ) > 
Joferyl ?-) is the s-wave spherical Bessel function and 
n is the principal quantum number. Due to factor- 
ization of the off-shell form factor in the half-off-shcll 
resonant amplitude, the ratio of the half-off-shell to 
on-shell astrophysical factors is S^ HOF ^/S^ ON ^ = 
k xA w (p xA ,k xA ( R) )/[(er 2v ^ A )Y xA (E xA )], where r\ xA is 
the Coulomb parameter. Taking into account that 
T xA (E xA ) - k xA P (E xA ), where P a (E xA ) is the barrier 
penetrability, we conclude that the energy dependence of 



the half-off-shell and on-shell astrophysical factors at low 
energies is the same. Note that if the Coulomb interac- 
tion in the initial state is included, the half-off-shell res- 
onant reaction amplitude should be divided by the same 
factor R as for direct processes. 

A simple model has been used to compare the energy 
behavior of the direct and resonant half-off-shell and on- 
shell astrophysical factors. The intention was not to re- 
produce experimental data, which requires a more so- 
phisticated approach, but to demonstrate that the en- 
ergy dependence of the half-off-shcll and on-shell astro- 
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FIG. 5: S(E) for 6 Li(p,a) 3 He S(E p6hi ) from the TH reac- 
tion 6 Li(d, an) 3 He (solid circles) [To| compared to direct data 
(open triangles [TH and open circles 0). The line shows the 
result of a second order polynomial fit to the TH reaction 
data |Jj}. 



physical factors are nearly identical when analyzed in the 
same model. Validating this makes it clear why the TH 
method is such a powerful indirect technique for nuclear 
astrophysics. The power of the TH method is seen in 
Figs 0] an d El ln both cases the energy dependence of 
the astrophysical factor determined from the TH reac- 
tion nicely reproduces the energy dependence of the as- 
trophysical factor obtained from direct measurements at 
higher energies. At lower energies the TH method pro- 
vides the astrophysical factor between bare nuclei while 
the direct data are distorted by electron screening. 
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